Entanglement of Trapped-Ion Clock States 
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A M0lmer-S0rensen entangling gate is realized for pairs of trapped ^^^Cd"*" ions using magnetic- 
field insensitive "clock" states and an implementation offering reduced sensitivity to optical phase 
drifts. The gate is used to generate the complete set of four entangled states, which are reconstructed 
and evaluated with quantum-state tomography. An average target-state fidelity of 0.79 is achieved, 
limited by available laser power and technical noise. The tomographic reconstruction of entangled 
states demonstrates universal quantum control of two ion-qubits, which through multiplexing can 
provide a route to scalable architectures for trapped-ion quantum computing. 

PACS numbers: 03.67.Mn,03.65.Wj,03.65.Ud,32.80.Pj 



Entangled states such as the famous EPR-Bohm states 
P, 0] have long been of interest in the interpretation of 
quantum mechanics 0; however, their generation has 
become a rapidly growing field with the recognition of 
entanglement as apowerful resource for quantum infor- 
mation processing!^. Laser-addressed trapped ions with 
qubits embedded in long-lived internal hyperfine levels 
hold significant advantages for quantum information ap- 
plications 1^. A critical issue is the robust genera- 
tion of scalable entanglement. In the context of trapped 
ions this is reduced to the problem of two-qubit entangle- 
ment, as plausible multiplexing schemes have been pro- 
posed to create a scalable architecture for a quantum 
processor^, 

Trapped-ion entangling gates mediated by phonons 
of the collective ion motion are susceptible to various 
forms of noise - qubit and motional decoherence, im- 
pure initial conditions, and technical issues associated 
with the optical Raman lasers driving the gateQ. Ro- 
bust schemes for gates based on spin-dependent forces 
have been proposed 0, 0, Q, 0| and experimentally 
implemented Il4| that, for example, relax the pu- 
rity requirement on the initial motional state of the ions. 
Here we report the realization of one such entangling 
gate for pairs of trapped ^^^Cd"*" ions that uses an advan- 
tageous implementationfisl of the M0lmer-S0rensen 
(MS) schemeHIIi. The implementation reduces sensi- 
tivity to optical phase drifts through an appropriate Ra- 
man beam setup and reduces sensitivity to magnetic field 
fluctuations through the use of magnetic-field insensitive 
clock states 0. 

Quantum state tomography 0, 0, |^ is used to 
characterize the gate performance for the creation of all 
four entangled Bell-like states. Previous applications of 
quantum state tomography with ions include the recon- 
struction of non-classical states of motion [2^ li^ as 
well as entangled states of optical ion-qubits composed of 
electronic levels |2^. Here we present the first such im- 
plementation for hyperfine qubits, in the process demon- 
strating universal quantum control of two clock-state ion- 
qubits. 



The MS gate for two trapped ions is based on opti- 
cal Raman couplings to the first vibrational sidebands 
of the ions' collective motion, assumed along the z-axis. 
The ions are equally illuminated by a bichromatic Raman 
field inducing simultaneous red and blue sideband inter- 
actions 22] that couple each ion's spin to the vibrational 
levels {\n)} of a single collective mode of motion - stretch 
or center of mass [Fig. ^a)] . With the Raman fields far 
detuned from the sidebands, negligible direct coupling 
occurs; however, the fields can combine to provide a res- 
onant two-step coupling between for example ||J,) and 
ItT), generating the entangled state l/%/2(||J,) + j|tT)) 
[Fig.^a)]. In general, the bichromatic field provides an 
entangling gate based on a nonlinear two-qubit interac- 
tion such as H ^ ~hfl{ax ^ ax)/'^, written in terms of 
Pauli operators. The coupling strength fi, which in the 
Lamb-Dicke limit is independent of the initial value of 
|n), is given by (r/f2)^/(5 where fl is the carrier Rabi fre- 
quency and S the detuning. The Lamb-Dicke parameter 
ry = kzo for the motional mode of interest is character- 
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FIG. 1: Two views of the M0lmer-S0rensen 'S' (Jx entan- 
gling gate for two ions in (a) energy space |3 and (b) motional 
phase space [3 for the gate-diagonal spin basis. The quan- 
tum number n and phase-space co-ordinates describe a given 
collective motional mode. Red and blue Raman sideband cou- 
plings are labeled by r and h and have detuning ih = &= —5r- 
For a closed phase-space trajectory, the phase "1> depends only 
on the area enclosed. 
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ized by Raman wavevector difference k along the z-axis 
of motion and zero-point wavepacket size Zo = y/h/2Muj, 
where co and M are the frequency and total mass of the 
excitation respectively. 

Reducing the detuning 6 accelerates the gate speed at 
the expense of populating intermediate motional states 
|25|. In this situation it is more natural to view the 
bichromatic field as generating a spin-dependent force 
constructed from red and blue sideband couplings with 
balanced Rabi frequencies and detunings 0, llfl. Efij |. 
The resulting interaction on each ion is equivalent to 
H ^ a^zFoSuiLOdt describing a (Jx-dependent force near 
resonance (a;d = w-f ^) with strength FoZo=r]^l. The total 
Hamiltonian is the sum of interactions on each ion. For a 
force resonant with the stretch mode and acting in-phase 
on the two ions, the time-evolution operator can be ex- 
pressed in the (Ta;-diagonal gate basis as a spin-dependent 
displacement as follows: 

= IT.T.XT.T.I + IUUXia.l 

where I?(a) is the displacement operator in the phase 
space of the driven normal mode [Figure 1(b)]. The 
value of the displacement is a{t,d) — ao(l — e^*"^*) and 
the corresponding phase accumulated over the trajec- 
tory is ^{t,S) — a^St — sinSt) in terms of the param- 
eter ao='r]fl/S. In general, the spin-dependent displace- 
ment entangles the spin and motional degrees of free- 
dom; however, for a closed trajectory {St = 27rm, m an 
integer), the spin and motion disentangle leaving only a 
spin-dependent geometric phase $g = 27rm(r/f2/J)^ ap- 
plied to the gate basis. A maximally entangling phase 
gate is constructed from a geometric phase of tt/2. We 
achieve this in the fastest time possible with m—1 requir- 
ing detuning S = 2r]n and gate time Tg = 27r/(5. Expressed 
in the computational basis, the MS gate structure, al- 
though slightly more complicated, makes the entangling 
action manifest: 

ltt>^*i = 7f (ltT)+*e'^=Ui)) 
|U)-*2 = ^(Ui)+*e-*HTT)) 

in)-*3 = ;^(in) + *e»*°UT)) 

UT)^*4 = ;^(UT)+»e-*^°|Ti)) 

The phases (f>o and (f>e have been included in the even and 
odd parity states, \E'i,2 and '^3^4 respectively, to account 
for the effect of both ac Stark shifts and Raman laser 
coherences, the latter modifying the spin dependence of 
the gate 0, For (j)^ = (j)^ — 0, the gate's action 

reduces to that oi a. Gx ® d'x coupling. 

Our qubit resides in the hyperfine clock states 
||) = |F=0,mp'=0) and ||) = |F=1, toj.=0) of a "iCd+ 
ion with frequency separation CL;h//27r=14.53GIIz and 
second-order Zeeman shift 600Hz / near zero magnetic 
field. A pair of ions, confined in a three-layer linear Paul 



trap [23 is aligned along the weak z-axis. The ions' sec- 
ular harmonic motion in the 2-direction is characterized 
by center-of-mass and stretch normal modes with fre- 
quencies Wc=2.05MIIz and LOg—^/^LOc- AH requirements 
for arbitrary two-qubit control are implemented as fol- 
lows. Both modes of motion are initialized to near their 
ground state [fie ~ 0.4 , ~ 0.2) with 60 pulses of Ra- 
man sideband cooling |23|. Due to the simplicity of the 
hyperfine structure (nuclear spin / = 1/2), the qubits 
are directly initialized with optical pumping to |t). Fol- 
lowing coherent operations the qubits can be read out 
with high fidelity using a photomultiplier tube (PMT) 
allowing |tT), and {|Ti), lit)} to be distinguished by 
virtue of state-dependent fluorescence Q. Unambiguous 
readout of all four qubit states, in particular and 
||t)j is achieved with an intensified CCD camera that 
can independently and simultaneously image the fluores- 
cence collected from each ion-qubit. A camera detection 
fidelity of 97% is achieved in 15ms limited by readout 
electronics. 

Coherent single qubit operations are achieved through 
a combination of applied microwave fields and ion- 
selective ac Stark shifts. Resonant microwaves provide 
simultaneous Rabi flopping of both qubits with a Rabi 
frequency of 56kIIz. Arbitrary independent qubit rota- 
tions are achieved by combining microwave operations 
with pulses of an off-resonant laser beam 200GIIz de- 
tuned from resonance. The laser beam with moderate 
waist (<10^m compared with the 2.5yLtm ion spacing) is 
aligned to be off-center with respect to the two ions, giv- 
ing rise to an intensity gradient and differential ac Stark 
shift between the two ions. A 10/iS exposure results in a 
phase shift difference between the two qubits of approx- 
imately 7r/2. 

The MS entangling gate is achieved using a pair of 
Raman laser beams 200GIIz detuned from optical reso- 
nance. An electrooptic modulator (EOM) at microwave 
frequency together with an acousto-optic modulator in 
each Raman beam provide the required bichromatic Ra- 
man beatnote [i^. The collective stretch mode is cho- 
sen for gate implementation due to the significantly sup- 
pressed heating rate [s^ . The stretch sideband Rabi fre- 
quency is typically 6kIIz. The wavevectors of the red and 
blue Raman fields are arranged in a counter-propagating 
geometry so that the spin coherence of the MS gate, in- 
cluded in 0e, is insensitive to optical phase drifts between 
the Raman beams [ll[3. Although this setup requires 
an accurate ion spacing to maximize gate speed, the sta- 
bility of the spin coherence is crucial for keeping the MS 
gate synchronized with the microwave fields during to- 
mography. A noise eater, stabilizing the Raman beam 
power, is used to suppress the effect of fluctuating ac 
Stark shifts (where the average qubit shift is 75kIIz com- 
pared with 13kHz gate speed as discussed below). 

The sequence to implement the gate begins with bal- 
ancing the sideband strengths to better than 10% and the 
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FIG. 2: Average brightness Sav (see text) versus MS gate de- 
tuning 5. Applied gate time (75^s) is within 10% of the ideal. 
Dotted line indicates expected signal modified to include an 
initial temperature fts =0.3^^. Solid line is a fit including 
offset and contrast factors to account for imperfections such 
as spontaneous emission. The fit gives a sideband Rabi fre- 
quency 77S7/27r=6.3kHz and initial stretch mode temperature 
ri,s=0.3. Vertical line shows ideal gate operation point S=2rifl, 
roughly at Sav = 1. Each point is the average of 150 PMT 
measurements. 
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FIG. 3: Time scan of MS gate showing (a) average brightess 
Sav and (b) parity 11. Ideal gate evolution shown as dot- 
ted lines with best fit including exponential damping shown 
with solid line. The fit gives a sideband Rabi frequency 
r;r2/27F=6.6kHz and detuning 5/27r = 12.8kHz « 2riQ/2n with 
other parameters the same as in figure |21 Vertical line shows 
gate operation time t = 2n /5~80fis. 

detunings to '--^lOOHz. Applying the bichromatic field for 
time T to the initial state |tT) while scanning the detun- 
ing S pinpoints the required gate detuning to near 2riQ 
[Fig. El . The dynamics of the frequency scan can be un- 
derstood in terms of the evolution of entangled states 
of spin and motion Assuming the initial spin 

state lit) and motional ground state \ns = 0), the aver- 
age ion brightness defined as Sav = 2P|| -I- P|| -I- P|| 
is SaviT,5) = l/2(l4-cos$(r,5)e-l"(^''^)l'/2). With the 
detuning now fixed, the average brightness is monitored 
while scanning the gate time. The time evolution reveals 
the overall spin dynamics modulated by faster dynam- 
ics associated with the phase-space evolution [Fig. |3I a)]. 
Each location of zero slope corresponds to the ion mo- 



FIG. 4: Parity versus phase of analysis 7r/2 pulse applied to 
the ^'i state. The solid line is a sinusoidal fit yielding an 
amplitude 0.79(2). The fidelity of the state shown is 0.83(2). 
Each point is an average over 50 PMT measurements and 
other parameters are as in text. 



tion returning on itself to form a closed trajectory. The 
return points are most clearly visualized in the parity sig- 
nal, n = (P^^+P^^)- (Pn+Pn) = 1/2(1 + e-2|"(-,«)r ) 
[Fig. ISIb)]. At the gate operation time (80/is), corre- 
sponding to the first return, the initial state |||) has 
evolved ideally to *i = 1/a/2(| TT) +ie*'^*=l ii))- 

The simplest indicator for the quality of the entan- 
gled states formed is the fidelity F— {^>\p\'i') with which 
the actual density matrix p matches the target state ^. 
The fidelity for creating the Bell-like states of Eqn. |21 is 
simply the sum of the two relevant diagonal population 
terms of p and the corresponding pair of off-diagonal co- 
herences. It is easy to directly extract the coherences for 
the even parity states ^'1^2 without single-qubit opera- 
tions. A single global tt/2 analysis pulse is applied to the 
state; varying the phase of the analysis pulse yields an 
oscillating parity signal [Fig. ^ with amplitude equal to 
twice the off-diagonal coherence 0, |^ . For the case of 
^1 as shown in Fig.^ a typical fidelity of 0.80 is achieved 
(which must exceed 0.5 to achieve entanglement [iMIs^ ). 

A full evaluation of the entangled state including a 
quantitative measure of the entanglement requires ac- 
cess to the full density matrix, in particular all the 
off-diagonal coherences. To determine the fifteen free 
parameters for a normalized two-qubit density matrix 
requires at least as many independent measurements. 
We follow closely the tomographic approach outlined in 
refs. 0, The density matrix can be decomposed 

in terms of a tensor product basis p = J2i j=o ^'ij'^i ® 
where ctq = I,ai = ax,cr2 = <Jy and 0-3 = Uz are the usual 
single-qubit Pauli matrices satisfying Tr{aiaj) = 2dij, 
and rij ~Tr{pai (g) aj) are real numbers. In the exper- 
iment we choose to perform projective measurements in 
the nine basis combinations {ai (g) aj, i,j = x, y, z} each 
yielding four possible outcomes for a total of twenty-seven 
independent measurements accounting for normalization. 
The fluorescence measurement accesses CTz projections. 
To implement transverse a^.y projections, we make use 
of independent single-qubit rotations to transform into 
the az basis before measurement. Repeated preparation 
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of a target state followed by tomographic measurement 
is performed for 200 shots per measurement basis. The 
total reconstruction time takes about 60s, dominated by 
the cooling cycle and camera readout time. 

A fast, direct inversion for the density matrix can be 
made with a minimum complete measurement set of fif- 
teen values Tij. However, this process in general leads 
to an unphysical density matrix due to experimental er- 
ror. Instead, maximum likelihood estimation is used to 
fit the data to a density matrix form constrained to be 
Hermitian, normalized and positive semidefinite. The 
inclusive and mutually exclusive nature of the four mea- 
surement outcomes for each basis is taken into account by 
least-squares weighting according to a multinomial distri- 
bution [s^ . Systematics of the tomographic process are 
assessed after the fact based on tomographic control runs 
of input states |tT) and assumed to be ideal. The 
results from the controls are used to extract detection 
biases (on the order of a few percent), microwave Rabi 
frequency and applied ac Stark shifts used for qubit ro- 
tations. Statistical errors for parameters calculated from 
the reconstructed density matrix are difficult to extract 
directly and so are obtained using a simple numerical 
bootstrap method The raw shot-by-shot data are 

randomly resampled with replacement to generate suc- 
cessive data sets from which a distribution of a parame- 
ter's value can be obtained. 

All four Bell-like entangled states are created according 
to Eqn.Elby applying the MS gate to the different compu- 
tational states. Figure |S1 shows their reconstructed den- 
sity matrices. The inferred fidelities for the target states 
^-i through ^-4 are F = {0.82(3),0.89(3),0.78(3),0.66(3)} 
where the phases 0e = — l.lrad and (/iq = 0.43rad are con- 
sidered free parameters obtained from the fits. The to- 
mographically obtained fidelity for agrees well with a 
simple parity-based assessment like that discussed above. 
The fidelity for creating the odd-parity states ^3^4 is 
worse because of inaccurate preparation of the input 
states Iti) and |jt) (i^~0.85). Accounting for this fac- 
tor, the fidelities of all states are on par. 

Inseparability (entanglement) of the reconstructed 
two-qubit states can be tested by performing a partial 
transpose of the density matrix and searching for a neg- 
ative value in the resultant eigenvalue spectrum [35L l3(il| . 
For example, the eigenvalue spectrum obtained for the 
*2 case is {-0.42(3),0.40(3),0.49(2),0.53(2)} compared 
with the ideal case {-0.5,0.5,0.5,0.5}. The negativity N 
[stI l3^ , twice the absolute value of the negative eigen- 
value, is obtained for all four target states V^i with val- 
ues iV = {0.74(6),0.84(7),0.60(5),0.42(6)}. Ranging from 
zero for a separable state to one for a maximally en- 
tangled one, the value gives an indication of the degree 
of the entanglement. Several quantitative measures of 
entanglement exist in the literature [13 j although lack- 
ing a closed form they are in general difficult to calcu- 
late. One standard measure that is directly calculable 
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FIG. 5: Tomographically reconstructed density matrices (a)- 
(d) for the four Bell-like entangled states ^'i through \l/4 as 
per Eqn. |21 To allow direct comparison of diagonal and off- 
diagonal elements, the reconstructed matrices were rotated 
into the real co-ordinate using fit parameter 0e = — l.lrad for 
(a) and (b), and 0o = 0.43rad for (c) and (d). Each state 
reconstruction uses 27 independent projective camera mea- 
surements averaged over 200 runs. 



for two qubits is the entanglement of formation Ep |40j 
again ranging from zero for a separable state to one for a 
maximally entangled one. In the context of pure states, 
the value of Ep can be interpreted as the number nEp 
of maximally entangled states required to reconstruct n 
copies of a given state |32j . The experimental values 
for the four states shown in Fig. [S] are Ep — {0.65(8), 
0.77(9),0.49(6),0.32(6)}. The entanglement of formation 
is a manifestly more strict indicator for the quality of an 
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entangled state than the fidehty and drops quickly with 
decreasing fidelity. 

Among the experimental sources of gate imperfection, 
spontaneous emission and fluctuating ac Stark shifts 
stand out as the likely primary sources of the observed 
infidelity. For our setup, the gate speed {Qg = ^ir/Tg = 
27] fl), which is proportional to the stimulated Raman 
Rabi frequency fi, scales as in terms of the op- 

tical linewidth 7/27r = 60MHz, the Raman laser intensity 
/ and its detuning A (not to be confused with the de- 
tuning d of the net two-photon Raman transition). In 
addition to generating the desired gate action, the Ra- 
man beams are responsible for a spontaneous scattering 
rate 7sc per ion and residual differential ac Stark shift Svst 
of the hyperfine qubit levels. The probability of a spon- 
taneous photon being scattered during a gate operation 
is Psc = '2"fscTg = 2Pj/ A where the factor of two accounts 
for the presence of two ions. A rough theoretical esti- 
mate can be made for the prefactor /? = •\/27''/eC'?~ 400, 
which includes a factor accounting for the bichro- 
matic field, a factor e ^ 0.2 characterizing the EOM Ra- 
man transition efficiency [i^, a Clebsch-Gordon related 
factor C = 0.5 and the Lamb-Dicke parameter 77 = 0.1 for 
the stretch mode. Similarly, the ac Stark phase acquired 
during a gate is (pst — SvstTg — f3ujhf / A. The relatively 
large value of u;?i//27r= 14.5GHz for Cd"'', while useful for 
high fidelity qubit detection, requires a significant detun- 
ing A to suppress Stark shifts. Experimentally, for our 
modest detuning A/27r«200GHz we measure a value of 
dvst/^TT = 75kHz, from which we obtain t/)^ — 127r and 
infer psc ~ 0.3. The value of Psc agrees roughly with 
the direct theoretical estimate and predicts an infidelity 
1 — F « 0.73psc = 0.2, roughly in agreement with the ob- 
served value for creating ^'1,2. The factor of 0.73 appears 
in the infidelity since a spontaneous scattering event will 
result in a mixed state that still has some residual overlap 
with the entangled target state. 

Increasing the detuning A can reduce the relative effect 
of both spontaneous emission and Stark shift (see also 
ref. however, a concomitant increase in the power of 
the Raman beams is required to maintain the speed of the 
entangling gate (and Raman cooling), thereby avoiding 
slower sources of noise such as magnetic field drift or 
laser beam-steering noise. In the short term, a reasonable 
increase in Raman laser detuning and power (currently 
~lmW) by a factor of ten would reduce the spontaneous 
emission and sensitivity to ac Stark shifts by the same 
amount. Ultimately detunings on the order of the large 
fine structure (74TIIz) of Cd+jl^ allow for significant 
suppression of both effects (see also ref. 0). 

In conclusion, a M0lmer-S0rensen gate has been re- 
alized to generate pair-wise entanglement of clock-state 
ion-qubits with reduced sensitivity to interferometric 
phase fluctuations of the Raman beams. The tomo- 
graphic reconstruction used to assess the resultant en- 
tangled states demonstrates universal two-qubit control. 



which is being directly applied to investigate prototype 
quantum algorithms |43|. 
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